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ABSTRACT 


We derive the closed form expression for the bit error probability of dense WDM 
systems employing an external OOK modulator. Our model is based upon a close 
approximation of the optical Fabry-Perot filter in the receiver as a single-pole RC filter for 
signals that are bandlimited to a frequency band approximately equal to one sixtieth of the 
Fabry-Perot filter's free spectral range. Our model can handle bit rates up to 2.5 Gb/s for 
a free spectral range of 3800 GHz and up to 5 Gb/s when the power penalty is 1 dB or 


less. 
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I. INTRODUCTION 


Wavelength division multiplexing (WDM) systems have been increasingly proposed 
as an attractive alternative to coherent optical frequency division multiplexing (FDM) 
systems [Ref. 1-5]. Although WDM systems with direct detection do not have the 
channel capacity of coherent optical FDM systems, they are much less costly to 
implement. Furthermore, present filter technology enables the designers to tightly pack 
the channel, resulting in dense WDM systems that can provide aggregate bit rates of many 
terabits per second (1 Tb/s = 10'* b/s). Dense WDM systems are particularly attractive in 
the area of undersea surveillance where hundreds of sensors and data collection sites are 
envisioned being merged onto single-fiber superhighways through massive data fusion. 
Other applications call for relatively low-bandwidth data collection over many months to 
be dumped quickly to a remote recording site in a matter of minutes. This 
“collection-and-dumped" compression can demand total data rates on the order of 
hundreds of Gb/s. Long distance between data collection sites and a remote recording site 
requires the use of optical amplifiers. Therefore, it is necessary to pack all channels within 
the optical amplifier bandwidth. 

A dense WDM receiver with on-off-keying (OOK) modulation can be modeled as 
shown in Fig. 1. Conceptually, the analysis involves two main operaticns: 1) a 


convolution operation to evaluate the signal at the output of the optical filter, a 








Fabry-Perot (FP) filter in our investigation, and 2) the integration of the output of the 
photodetector. Evaluation of bit error probability by the numerical analysis of these two 
operations has been carried out in [Ref. 6], with a number of approximations made to 
reduce the computational complexity. In this investigation the FP filter is shown to be 
well approximated by an RC filter within the frequency range | /- f, | < FSR / 20x, where 
FSR is the free spectral range of the FP filter [Ref. 7] and /, is the FP filter center 
frequency. For example, given FSR = 3800 GHz, the approximation works very well for 

| ff, | < 60.5 GHz; that is, the effects of adjacent channels within a 121 GHz bandwidth 
centered at f, must be included, while all others can be neglected. This simple model 
agrees well with [Ref. 6] as demonstrated in Section III. Furthermore, this model enables 
us to obtain a closed form analytical expression for the bit error probability for which 
numerical results can be obtained with little effort. Our investigation shows that this 
simple model provides accurate results as compared to those in (Ref. 6] for bit rates up to 
2.5 Gb/s when the effects of four adjacent channels are included with FSR = 3800 GHz. 
Actually, when the power penalty relative to single channel operation is 1 dB or less, there 
is virtually no difference in the effect of four or two adjacent channels. Thus, for this 
power penalty criterion, this simple model can handle bit rates up to 5 Gb/s for a FP filter's 
FSR = 3800 GHz. 

In Section II the closed form expression for the decision variable, and consequently, 

the bit error probability assuming all channels are bit synchronous as in [Ref. 6] is derived. 


Section III presents the numerical results which include the bit error probability versus the 





signal-to-noise ratio as function of the FP filter bandwidth and channel spacing, and the 
power penalty (relative to single channel operation without filtering or with filtering but no 
intersymbol interference) versus the channel spacing as a function of the bandwidth. 


Finally, a summary of results appears in Section IV. 








IL ANALYSIS 


The receiver model for the dense WDM system is shown in Fig. 1. The desired 
signal is filtered by a Fabry-Perot (FP) filter that rejects adjacent channels. The 
photodetector is assumed to have a responsivity @ (A/W). The detected current is 
amplified by a low noise amplifier that contributes a postdetection thermal noise n(/) with 
spectral density N, (A°/Hz). The decision variable at the output of the integration is 
compared to a threshold a to determine whether a bit zero or bit one was present. 
A. INPUT SIGNAL 

For convenience, we designate channel 0 as the desired channel, and channel k as 
an adjacent channel where &k = -M/2, .... -1, 1, .... 42 with M an even integer. We 


consider the equivalent lowpass (complex envelope) data signal in channel 0 and channel & 


as follows: 
bolt) = = boi pr(t- iT) (1) 
0 . 
b.() = Rs by,1e°*' pr(t — IT) (2) 
where 


T: bit duration 
5, = {0,1}: bit in channel 0 in the time interval (i7, (i+1)7) 


b, ,={0,e}is the /” bit in channel & in the time interval (/7,(/+1)T7) 
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6, : a phase offset between channel & and channel 0 and is assumed to be uniformly 
distributed in (0, 27) radians 
@, | radian frequency spacing between channel & and channel 0 with @, = -o, 
The function p,{t - i7) is defined as 
. 1, iT<t<(i+1)T 
Pr(t—iT) =o. otherwise (3) 
In both (1) and (2), the non-negative integers L, and LZ represent the number of bits in 


channel 0 and k, respectively, that proceed the detected bits 5,,. The received dense 


WDM equivalent lowpass signal at the input of the FP filter is given by 


rt) = JP bo Pb) a) 


where P is the received optical power. 
B. FABRY-PEROT FILTERED OUTPUT SIGNAL 
The FP filter can be characterized by the following equivalent lowpass transfer 


function [Ref. 1,7] 
1-A~p 





l- 
HS) = arm bd I-p 





I-p 1-A-p 
H ——?____ » 
(f)= I-pcos(2%)+jp sin(2&) Ip 


where p is the power reflectivity, A is the power absorption loss (zero for an ideal FP 
filter) and FSR is the free spectral range. For | f: < FSR/20n and assuming A = 0, we can 


approximate H (f) as follow: 


l- 
HS ay rae 


( —P)+/ FR 1WG= FRR 


lf |< FSR/20n 


—~4 


ce 


FSR(1-p) 


5 (6b) 


The free spectral range FSR can be related to the full width at half maximum (FWHM) 


bandwidth B and the finess F of the FP filter as 


FSR= "22" = BF 


Thus if the signal is bandlimited to | f| < FSR/20x, we can truly approximate (5) with a 


single-pole RC filter with the following transfer function and impulse response 





ee Ee 


h(th=ce“ , t>0 (9) 


Figures 2a-b show the magnitude and phase (radians) of H( f) of the FP filter in (5) and 
its single-pole RC filter approximation given in (8) for p = 0.99, F = 312.6, B = 12.1 GHz 
and FSR = 3800 GHz. Note that as the frequency increases, the phases of the FP filter 
and the RC filter differ markedly, but the magnitudes of their transfer functions remain 
identical and attenuate rapidly. When | / | > FSR/20r, the magnitude of H( / ) is very 

small, and therefore, the effect of adjacent channel interference beyond this frequency 


range is negligible. Figure. 3 shows the normalized impulse response of both FP and 





single-pole RC filters. In summary, the above approximation is valid for dense WDM 
analysis when the filter finess F is large or equivalently the FWHM bandwidth B is small 
since the equivalent lowpass signal must be bandlimited to about | f | < FSR/20n. 

This approximation has been used in [Ref. 5] to study spectral efficiency of optical 
FDM/ASK systems, which involves the evaluation of the decision variable for worst-case 
analysis using the eye diagram technique. Since we are interested in the detected bit 5 ,, 


in the time interval ( 0,7 ), we consider the output filtered signal s(t), 0< ¢ < T given by 


S(t) = Sa(t) + Sisi(t) + Saci(t), O<t<T (10) 


where 
5, (4): desired signal 
S,.,(t ): intersymbol interference signal 
8 
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Figure 2: Spectral characteristics of the Fabry-Perot filter and the approximated 
single-pole RC lowpass filter 








__ FP filter 


**** * approximate RC filter 





0 0.05 0.1 0.15 0.2 0.25 
time (ns) 


Figure 3: Normalized impulse response of the Fabry-Perot filter and the approximated 
single-pole RC lowpass filter 











Sa(t) = VP bo.0 h(t—t)dt 


= /Pboo(l-e“), O<t<T (11) 
pr 
si)= VP % bos [ Mt-s)ch 
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= /Pe~ z bo (et — eT), 0<t<T (12) 
i=-Lo 
ir. 
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0 
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k=-M/2 ry 
k#0 


—e(rion) TY) + by g(eV”" —1)}, O<t<T (13) 
The FP filtered output s(/) is detected by the photodetector which produces a 


current of @ | s(#) |’ Amps. This current plus additive white postdetection thermal noise 
current from the amplifier is integrated by the integrator to obtain a decision variable for 


the threshold detector. 











C. DECISION VARIABLES 
The decision variable Y appearing at the integrator output consists of the signal 


component X and noise component V 


Y=X+N (14) 
where 
rei &\s(t)|2at (15) 
0 
N -| n(t) at (16) 


We note that N is a zero mean Gaussian random variable with variance N,7. Substituting 
(10)-(13) into (15) we obtain the signal component _X as a function of the two parameters 
cT and @,7, which represent the effect of intersymbol interference and adjacent channel 


interference, respectively. 
X= #PTb,.(1- 201 -e*) + El -e 7) 


a ati pa e~2cT) [ 3 bo, (eT - elT))2 
i=-Lo 


M2 =! b 
PT ~-2cT kl 
+ €—(l-e | > ———_ 
er ( ) k=-M/2 I=-L 1+jo,%T/cT 


(e(¢TH0rTM +1) = e(cT+joeT)!)|2 





bt.obn.0 
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k#0 mre0 
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D. BIT ERROR PROBABILITY 


For a detection threshold a and an ISI/ACI bit pattern b = {5,,, 5, }, §=-Ly ...-1; 
/=-L, ..., 0; k= -M/2, .... M/2, k# 0; the conditional bit error probability of the OOK 
signal represented by the Gaussian random variable Y in (14)-(16) is given by [Ref. 8] 


P.(b) = 10> “at +3005) as) 
where Q (X) is defined as 
Q(x) = == | fe? "dy (19) 


and X, and X, are the values of X in (17) for 5,, = 0 and 5,, = 1, respectively. The 


average bit error probability P_ is obtained by taking the expected value of P, (6) in (18) 


14 


ss a ee ha a steel -_= eee ee saad 
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over all bit patterns 56. The minimum bit error probability is obtained by optimizing over 


the threshold a. In summary, given a, we calculate the following expectations: 
Pe =E {Pe(6)} (20) 


Pe, min = mink {P.(6)} =min sary p(b) 


M(L+1)+L, 
are 2M patterns 


where 


p(b) = 4b oe, gf ated yn. dOun 


X\ (O44... Or)-2 
+]- — ee) dD _ .a®D 
J re 1 INoT dD _yyp...dD yn} 
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il. NUMERICAL RESULTS 


In this section we present the numerical results for a) bit error probability versus 
signal-to-noise ratio Z = @P,J/T/No as a function of the normalized channel spacing 
(normalized to the bit rate) 


J=S =(Af)T (22) 


where Af= @,/ 22k is the equal channel spacing in GHz, and b) power penalty versus 
normalized channel spacing / as a function of the FP filter parameter c7= (x/ /p )BT 
where B is the filter full width at half maximum bandwidth(FWHM). 


Figures 4-5 show the minimum bit error probability versus signal-to-noise ratio 

(Z ) for cT = 5 and 10, respectively, and that of a single channel (SC) operation without 
filtering or with filtering but without /S/. In Fig. 4 we observe that a large degradation 
occurs due to /S/ for cT = 5 which represents a narrowband filter. As the FP filter 
bandwidth is made larger as in Fig. 5 with cT = 10, the ISI is reduced but the AC/ 
increases. 

In our model, we are constrained to M = 4 for the case under consideration. We set 
a FP filter with free spectral range, FWHM B=121.6 GHz, FSR = 3800 GHz, finess F = 


FSR/B = 312.6, and c = 38.4 GHz, 1/7 = 2.56 Gb/s, then c7T = 15. 
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Figure 4: Probability of bit error versus signal-to-noise ratio as a function of normalized 
channel spacing for ¢7=5 
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Figure 5: Probability of bit error versus signal-to-noise ratio as a function of normalized 
channel spacing for c7=10 
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(i.e. channel spacing is 12 times the bit rate or 30.4 GHz). In this model the farthest 
adajcent channel for M/ = 4 is twice the channel spacing which is 60.8 GHz. This verifies 
the assumption | /-f, | < FSR/20n = 60.5 GHz, where /, is the FP filter center frequency. 
This result agrees well with that in [Ref. 6, Figs. 6,9, M/F = 0.4,a=0.2]. Thus we 
incorporate the degradation caused by the four nearest adjacent channels. We observe 
that for bit rates of 1 Gb/s or less, our model is valid up to / = 10, and very little 
difference is observed between M = 4 and M= 10. Also, we observe that there is little 
difference between M/ = 2 and M = 4 when / > =10 for bit rates up to 3 Gb/s. In all results 
we set L, = 2 and L = 0. 

Figure 6 also shows the power penalty for a dense WDM system relative to a single 
channel operation at the minimum bit error probability of 10°. This is the required 
additional signal power (dBW) for the dense WDM system to be able to operate at the 
10'* bit error probability achieved in the single channel system with a SNR=12dB. The 
dense WDM system is /S/-limited at 2.2 dB, 1 dB, 0.5 dB, and 0.4 dB in power penalty 
for cT = 5, 10, 15, and 20, respectively. It is seen that for a 2.3 dB power penalty, the 
normalized channel spacing can be as close as / = 6 (i.e., a channel spacing of six times the 
bit rate) for cT = 5. If the power penalty criterion is 1 dB, the normalized channel spacing 
is /=12 for cT = 10, 15, 20. We remark that although the exact transfer functions of the 
FP filter is used in (Ref. 6], a number of approximation have been made to obtain 
numerical results. The approximations are 1) the /S/ is obtained by modeling FP filter as 


a single-pole RC filter [Ref 6, Eqs. (4) and (36)}, 2) approximating the finite integration 
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Figure 6: Power penalty versus normalized channel spacing as a function of Fabry-Perot 
filter parameter cT 











with an infinite integration in the calculation of AC/ [Ref. 6, Eq. (15)], and 3) the beat 
interference is ignored. On the other hand, the /S/ and ACI in our investigation are 
obtained by modeling the FP filter as a single-pole RC filter, using finite integration and 
including the beat interference. Since the results in our investigation and in [Ref. 6] agree 
well, we conclude that approximations are quite valid. We also note that our resuits also 
agree well with the simulation carried out in (Ref. 1, Fig. 17]. 

The above numerical results shown in Figs. 4-6 are obtained with an optimized 
threshold setting. Figure 7 shows the power penalty for fixed threshold a = @P7/2 which 
is the same optimum threshold for single channel operation (midpoint between the 
received power for bit zero and bit one). It is seen that the performance of a dense WDM 
system is quite sensitive to a for a narrow band filter. An additional 1.8 d3 is observed 
for cT = 5 for ] > 8, and 0.5 dB for cT = 10 for > 12. Negligible degradation is observed 
for cT = 15, 20 for J > 16. 

Figures 8-9 show the power penalty versus normalized channel spacing as a function 
of FP filter parameter c7 for the worst-case analysis with optimal threshold and fixed 
threshold, respectively. The worst-case bit pattern is fixed to produce the minimum 4X, 
and maximum X, where X, and X, are the values to X in Appendix-A equation (9) with 


b, , = 1 and b,, = 0 respectively. 
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Figure 7: Power penalty versus normalized channel spacing as a function of Fabry-Perot 
filter parameter cT with a fixed threshold a=0.5 
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Figure 8: Worst-case power penalty versus normalized channel spacing as a function of 
Fabry-Perot filter parameter cT with optimal threshold a=(X,+X,)/2 
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Figure 9: Worst-case power penalty versus normalized channel spacing as a function of 
Fabry-Perot filter parameter c7 with fixed threshold a=0.5 
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We observe that the power penalty for the worst-case analysis is only slightly larger 
than that of the exact analysis for / > 10 shown in Fig. 6. Similarly the power penalty for 
the worst-case analysis with fixed threshold is only slightly larger than that of the exact 
analysis with fixed threshold for / > 10 shown is Fig. 7. The reason for this is that for 
large channel spacing (/ > 10), the AC/ effect is small, so the AC/ bit pattern has a small 
influence on the power penalty. 


Figure 10 shows the normalized optimal threshold for the exact analysis shown in 
Fig. 6. It is observed that a ~0.4 for / > 10. Note that the normalized optimal threshold 
for the single channel operation is a = 0.5. 
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normalized channel spacing as @ 
function of Fabry-Perot filter parameter cT 


Figure 10: The normalized optimal threshold versus 
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IV. CONCLUSIONS 


We have presented a simple model for the analysis of dense WDM systems 
employing an external OOK modulator. The only approximation that we use involves the 
modeling of the Fabry-Perot filter by a single-pole RC filter assuming the equivalent 
lowpass signal is bandlimited to the frequency range | /|< FSR / 20x. This model 
enables us to obtain a closed form expression for the bit error probability which previously 
can only be obtained via numerical analysis (Ref. 6]. For FP filter with an FSR around 
3800 GHz, our model can include the AC/ effects of four adjacent channels for bit rates 
up to 2.5 Gb/s. Our numerical results show that this model agrees well with that in 


[Ref. 6]. 
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DERIVATION OF FORMULA FOR DECISION VARIABLES 


DESIRED CHANNEL : CHANNEL 0 
ADJACENT CHANNEL : k, k = —M/2,...,-1, 1,...,M/2;M: even 


BIT IN CHANNEL 0 IN ith TIME INTERVAL (i7,(i+ 1)7) - 
bo; € {0,1}, boo : DETECTED BIT IN (0, 7) 
BIT IN CHANNEL k IN /th TIME INTERVAL (IT, (1+ 1)7) ; 


by; € {0,e**} WHERE j = J-1_ - IMAGINARY NUMBLER 


@,: FREQUENCY SPACING BETWEEN CHANNEL k AND CHANNEL 0, 
O, = -O, 

¢, : PHASE OFFSET BETWEEN CHANNELS k AND CHANNEL 0, 

UNIFORMLY DISTRIBUTED BETWEEN (0, 22) 


DATA SIGNAL IN CHANNEL 0 : 


bo(t) = % Bospr(t- i 
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DATA SIGNAL IN CHANNEL k : 


bi = by bye" pr(t- IT) 
L,,L: INTEGERS 


1, O<r<T 
Pp r(t) = { 0, otherwise 


wan 1 IT <t <i4+)T 
Pp r(t —i T)={ 0, otherwise 
THE RECEIVED SIGNAL AT THE INPUT OF THE FP FILTER OF 


CHANNEL OIS : 
rp) = JP bolo Pb 


P: RECEIVED OPTICAL POWER. 


THE OUTPUT OF THE FP FILTER IS 


oc 
ro(t) = ne h(t — t)r(x) at 
WHERE /(#) IS THE EQUIVALENT LOWPASS IMPULSE RESPONSE OF 


THE FP FILTER OF CHANNEL 0 


ro)= WP f ae-abooar+ JP Ef me—bu(r)er 
ee ke0 
= IP boo J h(t-t)pr(t)dt 


+/P = boi I A(t — t)pr(t — iT)dt 
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+/P ne > bes I h(t — re" p r(x — [Td 


SINCE THE DETECTION INTERVAL IS 0 <#< 7, WE ONLY NEED TO 


EVALUATE s(t) =7,(t), 0 <#<T 


t = (i+1)T 
s(t) = JP boo | h(t-t)dt+ VP . % boi J A(t-dat 
0 i=-Lo T 


i 


= (i4+1)T 
WP SY gd oy f ne—demdt] 
oe l=-L IT 


+bx0 h(t — te? dt} 
0<t<T 
=s,(t)+s,.At)+s ,-At) (1) 
S,(t) : DESIRED SIGNAL 
S,sAt) : INTERSYMBOL INTERFERENCE 


Sicht): ADJACENT CHANNEL INTERFERENCE 


FP FILTER : 
LOWPASS EQUIVALENT TRANSFER FUNCTION 
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1-A-p _ I-p 1-A-p 
Hf) = — tae op ee 
ee \-pcos( 2 jp sin( 2) 


p: POWER REFLECTIVITY 

A: POWER ABSORPTION LOSS (A= 0 FOR IDEAL FILTER) 
FSR: FREE SPECTRAL RANGE 

SINCE {<< FSR (FOR OPERATING FREQUENCY RANGE) 


WE CAN APPROXIMATE H(f) AS (ASSUME A= 0) : 





1 
H(f)= ares 3 
f (py 1- ve 1 ee 
FSR(1-p) 
Hf) * a9 re WHERE = 
FOR FP FILTER WE ALSO HAVE 
FSR _ "JP 
B~ I-p 


B. FULL WIDTH AT HALF MAXIMUM BANDWIDTH (FWHM) OR HALF 
POWER BANDWIDTH 
THUS 


ce“ ¢>0 
h() = Lo, otherwise (2) 


TION | _19 


1. S, (0): 
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O0<t<T 


f h(t -1)pr(x)dt 


= h(t — t)dt 
0 


h(-d) 





t>T 


f h(t-1)pr(at 


T 
=| hit = t)dt 
0 





t 
sa() = VP boo J h(t—dadt O0<t<T (3) 
0 
SUBSTITUTING (2) INTO (3) WE OBTAIN 
t t 
sa(t) = JP boo | ce dt = VP boo ce | eat 
0 0 


= JP booce “tet l) = JP booe“"(e* — 1) 
= JP boo(l-e~“) 0<t<T 
a) BIT 1: sp} (ft) = /P(1-e~“) 0<t<T 


b) BITO: spo(t) =0 0<t<T 
2. Sis, (f): 


pr(t—iT) 


‘Secdhng 


iT (i+ 1)T 





f A(t —t)pr(t — iT)dt =0 
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iT<t<(it+1)T 11> 4 DT 


SS 


iT t G@+1)T 





[ he-vpr(e- iT)dt [ ha-vpr(e- iT)dt 


t (i+1)T 
=| h(t—1)dt = | ht-dt 
iT iT 
= (i+1)T 
sis) = JP = bo, | ht-vdt O<t<T (5) 
i=-Lo iT 


SUBSTITUTING (2) INTO (5) WE OBTAIN 


sis(t) = JP & do [coat 


-} ae 
iT 


-} mare 
=-/P > boi cet 


i=-Lo 


sisi(t) = JP z bos ee" XT — e*7) 
i=-Lo 
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sis(t) = Pe = bo (eT -—e*T) O<t<T © 
i=-Lo 


WORST-CASE ISI: 69; =6- , Lo = 
Sis) = JP be 0<t<T (6a) 
3. Sacit) : 


pr(t—/T) 


ce 


IT (/+1)T 
/=-l,...,-1 





SO THIS IS THE COMBINATION OF THE ABOVE 2 CASES 


. : (+41) | 
sat) =vP "© ot by f AG—de™ ar] 
Pane l=-L IT 


t 
+byo | h(t—t)e!*dr} O<t<T (7) 
0 
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SUBSTITUTING (2) INTO (7) WE OBTAIN 


JP + 3 b ea —ft-t) pjoet fe 
P—J e = k 
Sact(t) = JP bien 2, kl , C ] 


t as 
+bxo | ce el dt} 
0 


(/+1)T 


mp =i | 
= JP ¢ {LL berce* fF etoutgy 
ter l=-L IT 


t 
+by0ce-* [ eld} 
0 


= JP e 129 byice 1 RH TY 


kM? C+jOr 


—cte elVOKe yy 


+b kO ce C+IO, o} 





2 
Saci(t) = JP ce~* 2 SLE a b; eCvoxXR DF _ elctarylT)] 
k+0 


+b, (eo! = 1)} O<t< T (8) 


WORSE-CASE ACI : 54) = bk = 6, L=0 
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2 
we beCven)! 
Sacf) = JP ce 2a rte 
k#0 





[2 
sich) = JP b a a 0<t<T (8a) 
k+0 


OUTPUT OF PHOTODETECTOR: @ |s()| ¢ 


T 
SIGNAL OUTPUT OF INTEGRATOR: X= #f |s(1)|7dt 
0 


\s(2)|? = Isp(t) + sisr(t) + sacr(? 
= |sa(|? + |sis(O|? +1sacr(Ol? + 2Re{sa(sjs)(0)} 
+2Re{s3(O)saci(t)} + 2Re{Sis(Dsact)} 
WORSE-CASE 
1) sa() = VP boo(1 -e“) 0<t<T 
Isa()l? = Pb3 (1 -e)? = Pho o(l-2e“ +e") O<t<T 


t 2 2 2 t 2 t 
J [sa(t)|"dt = PTbo 9 — 2Pboo J e“dt+ Pboo } edt 
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= PTb?  - 2Pbpo-%e~'l5 + Phoose "15 
= PTbjy + 2Pbi9 (9 -F -1)- Pio (eter -1) 
= PTb) [1 + (e*? - 1) - se? - 1)) 
= PTb} [1 - (1 -e~*7) + (1 -e 7) 


2) s¥o(t) = JP b-e~“ O0<t<T 


|s"5(0)|° = Pb2e-2" O<t<T 


t we 2 2 t 2¢ 2_1 2ct|T 
J Isrq(Ol? ¢ = Pb2 J edt = Pb2-—-e*""| 


= —PH(e-2eF _ 1) = FE] - e¥7) 


PTb?2 2 
=i or ( l-e 2cP) 


w 2 Ke 
3) sedd= PbS Se O<t<T 
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styl? = Pio? “eS EM ctr 
b=-MN2 MP2 (14738 1-752) 


T 2 
we 2H 2 —_i— (@ pO») 
J IsxcxO| Oe Pipl in pe (14st 1-72) . a 


2 
= Pj? “EP 4 tt oflore, 7 
| | k=-M/2h=-MI2 (14jEX1-j2)@ tn ) 
ke he0 


e 1 
as P\b b|? 


¢ 


rin heh ow why (=t-=A) 
[eter | 


pro? M2 Me seb 
k=-MI2 h=-MI2 W14yEX 1-2) 
k#0 = h#0 





@, = 2, 1= INTEGER > 0 


t (@4—-On T, ©,=@p 
THEN Je! k-Oh dt = {0 wis, 


[ [swe ¢o[2ae = Priel? —— 
0 


k=-MI2 1+(-4)? 
ke 
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4) 2Re{sp(t)siq(0)} = 2Re{ JP boo(1 —e“") VP b-e} 


= 2Pbo9 b-(e“' —e-*") O<t<T 
2Re{sa(t)sia(t)}dt = 2Pboob-{<=13 — 17} 


= 2Pbo0b-{<—— + —! 


PTbo.0b- 


= (1 - —2e~7 + @-%T) 


s) 2Re{s*(t)saci(t)} = 2Re{ VP boo(1 -e) VP b , 


Op 
k=-M/2 lt 
k#0 


/2 kt —p (o-p) 
=2PhopRe{b & See) 
Mn Ve 





l 
| aRe(si(osscx}de= 2PbhooRe{b ie ye 


T T 
[J edt-[ edi} 
0 0 


0; = a I= INTEGER>0 














T 
Ort de — ajo — 1) = 
Je’ dt = —~(e! 1)=0 


Biba aie l 11 ~cT gfe, L l-e? 


_l leew? 
1 ) 


2Re{s*(t)saci(t)}dt = 2PbooRe{b in 








= 2ProesRetthy _ e~T) +> 1 
k=-MI2 1458)? 
ke0 


6) 2Re{ste(i)Saci(t)} =2Re{ VPbe (Pb % Se 


kM 1+/F 
ke0 


a 


= 2Pb_Re{b aa TRE? 


| 2Retsiglsacx()}= 2Pb-Re{b} in eo 


= PRM -e~7) eae 


Mn 148)? 


4) 





SUMMARY OF WORST-CASE ANALYSIS 


@,= 72 — 1 =INTEGER>0 
T 2 

X= 2@ { \s(o|"dt 
0 


2 
= # {PTb? {1 - 2(1 —e*7) + tel —e 7) + FEI - 7) 


+PT|b|? - 1 + Peeh(1 —2e~T +e) 


ke-M2 np 
ke0 


2PTbooRe{5} (yo -cT } 
aie ai HF)? 


pre -cT 1 
re in wry 


X= #PT{b3 [1-21 -e*") + sl -e *)] 


+b2[=te(1 — e-%7)] + Boob-[-p(I - 2e*7 +e) 


+{|bl? + (0 —e-*7)Re{b}(b-- boo] & 


1 
k=-Mi2 142)? 
k#0 
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WHERE b € {0,e%}, Re{b} © {0,cos}, |b]? © {0,1} 


EXACT ANALYSIS 
T 
nf Isa(O|2dt = PTB? of! - 2 - e~7) + (1 - e-**7)) 
0 
2) sis(t) = Pe : bo,(e*? — a7) O<t<T 
Isin(|? = Pe] B boyle*T—eM)}? Oct T 
1=~L0 


t x ; 
J [sis()|7at = 21 ~My + ; bole"! — eT)? 


2 
3) Saci{t) = IP ce ¢ Le —{[ z by Ke vex Ler 


k=-M/2 Clty) l=-. 
k#0 


eli P)] + by (el YF" — 1} O<t<T 
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Saci{t) = JP e~! in 2 x - sp hadel et ter 


ell Eel) 4 JP JP ini = ater" - e~) 


0O<t<T 


~— by eve tM l+lcT _ ells eT |2 





Isaci()|? = Pe~™*| ¥ 


k=-M?2 Is-L l+j 
ke 


+P +> Me — ante = en (eJon! — e~) 
sr casei (4ySK 1 


x bib, Sk 
42PRe¢ “SOB etme _ (et xtrnver 
k=-M?2 m=-M?2 l=-L (14j-EX 1-72) 
k+0 m+#0 


St om rs 
—e lS VeTyeXtly =)ct -@ 2cty) 


r 
EVALUATE | |saci(#)|7d¢ : 
0 


ZI = 
a) j Isaci(t)| {at = mal — e~*T)| + : atl. 
0 


k=-M)2 I=-L 1 4+)-F 
ke0 


(ells EM eT =. ellaPylery| 2 
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Len 


SPECIAL CASE: @; = ou 1 = INTEGER > 0 


mR =w3 b 
= 20) _¢-2Ty| Mt 3 I_(p(l+l)eF _ gleTy|2 
_- —-@ e —-@ 
<r ) k=-M2 I=-L 14/2 ) 


b b° " 
7 2 Me os k0° mo | —@ ¥, d 
{ pio. im t 


T T . T 
—f ees gt— f etiam t+ | e-2 dt} 
0 0 0 


bob mo T, Q=Om 
Op 0 feK-om)T_ 
ee a (lysF) aan , 0, FW» 


~ TL jeom )T _ 
ete ba eoee + 2(1-e7)} 


C-JO, C+fO m 
biobs, = 
St A ive et 
cT Ato) ee 
Temes «2 m 


k= M2 mM (ayy) 


@ 
Kg) 


jeer 14 eT _ 
ets re Ka: ] +i(1-e*7)} 


1-4 ly? 


SPECIAL CASE: @; = 48 1 = INTEGER > 0 
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T eT 
F a = T, @c=Om — 5 ©: Ok Om 
frelerteldt = Geos = 45 ayaa, 
eT — g-jorT = |} 
=p Tt us _Lbeol?_ Me brobuo esol 
bean 2)? +e in m-MN2 (14j252 (1-28) * 1-288 
—~<T_y as 
: x(1 -e-**7))} 
bub 
Sy) 


——o—_ 
i 


C) Isac()|5 dt = = 2PRe{ in nr=—Mr2 I=-L (14+j 


(ev EXl+ler = e(lveWeD yy j ely yet yy edt} 
bk,10 mo 


12 
=22 Ret "§ 
b= M2 m=—M2 I=L (at y 1) 
m<#0 
(el vPXE Der a eller Ever ype a 1 - —e-XT)}} 
oa 1 = INTEGER > 0 


SPECIAL CASE: @, = 








Ss brid 
= 25 Ret + 6) z aap ee 


be MND m= MP2 Ie-L (1421-22) 
ke0 
-}(1-e*7)}} 


r -— 
| Isacr Oat = EEC - e-*7)I + SPL (ell exter 
0 


b=-Mi2I=-L Leja 


etter 4 2 PG eee 
k#0 m+#0 


o 
ce%, omen pg 
a 
oes yt 
oO, Om » OFOm 
Kore) 


om 
em )eT_] 


lyy _ 4-2cT 
espe tle 


2 <li, 
+27FRe{ 2 ev —— 
Fir hale fa le-L (14j-eX 1-74) 


om 
(e(luENlH eT - elluPueT yy j—e “eT 


ly 


“Hl -e*7)} 








SPECIAL CASE ° 


T = b 

24, _ PT -2cT | 3 kl (a (I+leT _ 4hcT \2 
f lsacr(O|"dt = (1 - e**") ¥ aw (e e™") 
5 2-7 k=-M?2 I=-L 14j2 


+PT{ Se toeol 2 Me brobino 
kez +2)? zt k=—-MI2 m=-Mi2. (1+)72\(1=/22 
ke0 ke ma0 





-cT_ e-T-] 1 = 
(a + er +3 -e*7)} 
cT cT 


wre FS 2 amer™ 


Midm=—MI2i=—L (142821 -28 
“ge; 


l-e7*7 


<eS 4/2 - (1 —e*T)}} 
T 
4) | 2Re{sa(Osist(4) dt 
ef 2[ VP boo(1 —e)] [VP e™ 2 bos(e"*'*" — e'*7) dt 
0 I=-LO 


= 2Pbho0 z boi (etter ~ ell) (et — ev dt 


i=-L 0 


Cove beet nein innit Rize eniafe cenntccerdaten at etelenistins cera seitetihn Sed gehts eel 











T 
J 2Re{sp(tSjs)(t) }at 
-) ; . ent = 
= 2Pboo 2 bo (eT — eT S15 -£ — =|0] 


= “hee Dy bo i(elrbeT — eT) - eT + Se Ur _ 3] 


j=-L 


_ PTboo 2 bo (eT — eT (1 +.e-%T —2e~T) 


i= 


a Freon) ae e-2cT c 2e~') 3 bo (eT - ell) 
cT i=-Lo 
r * 
s) | 2Re{sz(saci()}at 
0 


=2 [ Ret{vP bool -e~)) 
0 


2 b 2k 4 
[/Pe* me. ; te £\1+1)cT _ eilwz MeT) 
K20 


ee ge er(ere — e-*)] at 


ken V+) 














L 
=2 | Re{Pbhoo(1-e“)e™ 
0 


R <i (ly EM+DeT _ el vZVeTy + Pho o 


k=-M/2 l=-L He 
k#0 


2 be —te_(er -~ete el ~j£)et es e-2) at 


k=-M/2 i+jt A 
k+0 


5 2 Sk had 
= 2PbooRe{ . ¥ —(el4 PMT _ p (lyse) 
k=-MiU=-L +i 





r 12 
i} et — e 2") dt+ + mt 


r ° .Op 
j (elon! — et eae + edt} 
0 


2 = 
= 2PbooRe{ 23 + 5¢-%*" — e“7) FS Pi (g(yF MoT 
kMdaL We 


ee exT-1 eo *T-1 +e) 


( 1+j—£yleT 
~~ + & JOx ¢ 2c (1 ~j Le 


k=-M/2 ly 
k+0 
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2 - 
_ Poo Re (j+e%T eT) “PY Pe (elt Entever 
cT ben L. lee 


Ok 2 bk ek] oe 1 
—e(lewel) +2 ¢ (a +e -1 —35e ul +S 
sare lj IS 


SPECIAL CASE : ©; = a 1 = FIXED INTEGER > 0 





r . 
J elrtdt = <-(e7 - 1) =0 


THEN 
T 
! 2Re{sp(t)saci(t)}at 


= “T*Re{(1 + e@2eT _ 2e~T) Mf kl me 2 yl+1)cT 
k=- thas 1422 


oS 


a 


-2Rkl 
~e (rel) 4.2 ¢ bko —<(e~7 - [= teu 4 + +4 
k=-M/2 142 
ke0 











T 
J 2Re {$75)()Saci(t) at 


= 2Re | (NP Pe z bo (et? — e*7)] 


i=-L 


a b Ok bad 2 
(VP e~ i) ale e(lyz A)(41)eT _ p(lti7 McT) 
k#0 


4: VP 2 M2 bio ar(emt — e~')ldt 


= b ib ; 
= 2PRe{ _ = » yz 0, ror ails ell) 


2 — 
(ely EneI)e _ elltiPyeTy Je*ate Y ! bosbko 
k= MN t=-Le Lajt 


T r 
(e(tDeT _ eicT) [ (eter —e~*")dt} 
0 


S  boib 
= 2PRe{ 2 3 3 0 —(e (i+1)cT _ gic?) 
Mi2 1=-L iz-Lo 1+ 


(el eXl er = elle ery L(y _ e~2T) 
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2 =) bob iy Per = 
3 0. tabi (gieer _ ge peer +(e 2cT _ 1)}} 


ak 
b= M2 i=-Lo layst tye 
ty bo (e"*!*T — eT)Re if 3 1-—e-*T 
24, il )Re { olin EL yt ) 


why 
(ev EXI Der = e(laziel), Me bro je le +e 72eT _ 1)} 
kM lj lj 


SPECIAL CASE: @; = oat! ] = INTEGER > 0 


T 
J 2Re { Ss Isi(f)SACI (Adt 


ee 4 5 3 bo (i+1)cT icT 3 bk -2cT 
= e —e')Re ¢ l-e 
cT 5-7, al Rei k=-MI2 I=-L 142 Tea ) 


1 I+lyeT _ IcT bro ¢2-2e~ a 
(e(lEMtDeT _ Q(t ye oe wey (OES a) +e UT _ 4) 
k+0 


WHERE e(ltyAMl+ eT = e(ctio M47 — g(cT+joxT\i+1) 


aad eo T(+1) gjaeT(+1) = ecT(+1) 














SIMILARLY=> eV EVET = elcviasNT — gclT 


T 
=> 2Re{Sis(t)Sacr(t)} at 


EE bole? —eT)Rel SS ee 


T 4 eM Fa-L 14j2 





(1 — 2-27 (eeT+1) — eT) Me bua als) + ~4eP _ $)} 


k=- rr 1428 

















APPENDIX B 


MATLAB COMPUTER PROGRAMS 


exact analysis with optimal threshold signal out of the integrator, the signal X contains desired 
signal and ACI & ISI and postdetect noise , The formula for the BER is - 


pd) 


i 
Misia 2M pattern 


where p(b) = HgLyMf f -- ij | Oat dO un ...d®un 


Sie fo OC TS dO un. } 


M is number of adjacent channels. Since we assume upper and lower channels are synchronized 
individually, i.e. for our model the power of (1/27) is fixed to 2, also we only have two 
integrations and two arguments. 


M=4, k=[-M/2:-1 1:M/2]; 
m=k; 

% produce the controlled matrix b to control 64 different bit patterns 
m1=[ zeros( 1,32) ones(1,32) J; 
m2={ zeros(1,16) ones(1,16) zeros(1,16) ones(1,16)}; 
m3=(]; 

for i=1:4 

m3=[m3 [zeros(1,8) ones(1,8)]]; 
end 


m4=(], 
for i=1:8 
m4=[m4 [zeros(1,4) ones(1,4)]}; 
end 


m5=zeros{ 1,64); 
m6=(], 
for i=1:16 
m6=(m6 [zeros(1,2) ones(1,2)}]; 
end 
m7=(), 
for i=1:32 
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m7=(m7 zeros(1,1) ones(1,1)}; 
end 


b=[m1;m2;m3;m4;m5;m6,m7]', 

% signal to noise ratio range in dB 

RPTN_DB=(10:0.5:25]; %RPTN_dB 
ppp=10.(0.1*RPTN_DB); 

len!=length(ppp), 

% solalph function is provided by Professor Randy L. Borchardt 


n=10; 
[bpx,wfx}=grule(n), “%bpx=bpy ,wfk=wfy 


% single channel 
BERO=0.5*erfc(ppp/8“0.5); 


pp=(]: thresh=[]; thresh is not normalized threshold 


for CT=[5 10 15 20] %cT is Fabry-Perot filter parameter 


if CT==5 
qqq=(2:12 17 20}; 
elseif CT==10 
qqq=(3:12 17 20); 
elseif CT==15 
qqq=[4:12 17 20); 
elseif CT==20 


qqq=(5:12 17 20); 
end 


pe=[]; thresh1=(]; 


for I=qqq 
BER=(], thresh2=[), 


for RPTN=ppp 
ap=linspace(0,182,11); | Yeapproximated thresholds 
x3=0; % first few loops to find out the threshold which make the 
% BER minimun don't care about the scale 


for i=1:64 x3=x3+... 
solalph('xx00',0,2*pi,2,bpx, w£x,0,2*pi,2,bpx, wfx,ap,b,CT,i,1,k,m,RPTN)¢... 
solalph(‘xx1 1',0,2* pi,2,bpx, wfx,0,2* pi,2,bpx, wix,ap,b,CT, i, .k,m,RPTN), 
end 
{val,ind]=min(x3),; 

















lef=ap(ind)- 16.6, 
if lef<0 lef=0; end % to aviod the threshold go beyond the negative side 


ap=linspace(lef,ap(ind)+ 16.6, 11); 
% two more time to find alpha 
x3=0; 


for i=1:64 x3=x3+... 
solalph(‘xx00',0,2* pi,2,bpx, wfx,0,2*pi,2,bpx, wix,ap,b,CT,i,..k,m,RPTN)+... 
solalph(‘xx1 1',0,2 *pi,2,bpx, w£x,0,2*pi,2, bpx, wix,ap,b,CT,i,.k,m,RPTN), 
end 


(val,ind]=min(x3); 
lef=ap(ind)-3.1, 


if lef<0O ,lef=0; end 


ap=linspace(lef,ap(ind)+3.1,11); 
% three more time to find alpha 
x3=0; 


for i= 1:64 
x3=x3+... 
solalph('xx00',0,2*pi,2, bpx, wix,0,2*pi,2,bpx, wfx,ap,b,CT,i,I,k,m,RPTN)+... 
solalph(‘xx1 1',0,2*pi,2, bpx, w£x,0,2* pi,2,bpx, wfx,ap,b, CT, i,1,k,m,RPTN), 
end 


[val,ind]=min(x3); 
lef=ap(ind)-0.57; 


if lef<O ,lef=0; end 


ap=linspace(lef,ap(ind)+0.57,8); 
% four more time to find alpha 
x3=0, 


for 1=1:64 
x3=x3+... 
solalph(‘xx00',0,2* pi,2,bpx, wfx,0,2* pi,2,bpx, wix,ap,b,CT,i,.,.k,m,RPTN)¢... 
solalph(‘xx1 1',0,2*pi,2,bpx, wfx,0,2*pi,2, bpx, wfx,ap,b,CT,i,I,k,m,RPTN), 
end 


[val ind}=min(x3), 
ap=ap(ind), 
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% after find optimal alpha use double integration 
qq3=0; 


for i=1:64 
qq3=qq3+( dbgquadm('xx00',0,2*pi,2,bpx, wix,0,2*pi,2, ... 
bpx,wfx,ap,b,CT,i,..k,m,RPTN)*... 
dbgquadm(‘xx1 1',0,2*pi,2,bpx, w£x,0,2°pi,2..... 
bpx,wéx,ap,b,CT,i,.k,m,RPTN) )/(1024*pi*2), 
end 


BER=(BER qq3]; thresh2=(thresh2 ap]; 


if qq3<10“(-15) %for save time only interesting in 10“(-15) 
ii=find(ppp==RPTN), 
BER(ii+1:len1)=5* 10(-116)*ones(_1,length(ii+1:len1) ); 
thresh2(ii+ 1 :len1)=(ap+2)*ones(1,length(ii+1:len1) ); 
break 
end 
end 


pe=[pe;BER], thresh1=(thresh1 ;thresh2)],; 
end 


See lh thresh=[thresh;thresh 1], 


time2=toc; 
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% worse case form appendix setting optimal threshold equal (x0+x1)/2 

M=4,; 

I=linspace(0,6, 101); % set I value in linspceinteger to find minimum I 
RPTN_DB=0:0.2:20, % x-axis dB range 
RPTN=10.(0.1*RPTN_DB);, % change to ratio 

% single channel 

BERO=0.5*erfc(RPTN/8“0.5); 

Yefind optimal alpha 


p=(]: 
for CT=[5 10 15 20] 
x3=0, 
for k=1:0.5*M 
X3=x3+2./(1+(2*pi*k*V/CT).%2); 

end 
x0=(1-exp(-2*CT))/(2*CT)+x3*(1+2*(1-exp(-CT)VCT); 
xl=1-2*(1-exp(-CT))/CT+(1-exp(-2*CT))/(2*CT); 
[val,ind}=min(abs(x0-x1)); 
I=ceil(I(ind)); 

I=(I(ind)); % our minimum I value 


% set different I and to find the BER 
ss=]+1-I+10; 
BER=(]; 
for I=ss 
x3=0; 
for k=1:0.5*M 
x3=x3+2/(1+(2*pi*k* /CT)*2); 
end 
x0=(1-exp(-2*CT))/(2*CT)+x3*(1+2*(1-exp(-CT))/CT); 
x1=1-2*(1-exp(-CT)VCT+(1 -exp(-2*CT))/(2*CT), 
BER=[BER;0.5*erfc( RPTN*(x1-x0)/8%0.5 ) J; 
end 
pp=[pp:BER]; 
end 


semilogy(RPTN_DB,BERO’'--,RPTN_DB,pp(:,:)) 


axis([10 19 10°(-15) 1]) 
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